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THE SADDLE-POINT METHOD AND THE 
LI COEFFICIENTS 

KAMEL MAZHOUDA 


Abstract. In this paper, we apply the saddle-point method in conjunction with the 
theory of the Ndrlund-Rice integrals to derive a precise asymptotic formula for the 
generalized Li coefficients established by Omar and Mazhouda. Actually, for any 
function F in the Selberg class S and under the Generalized Riemann Hypothesis, 
we have 


Air(n) 



log n + cpn + 0 (\/n logn), 


with 

cf = ^(7 - 1) + bog(AQ|,), A = ]^ A^^^ 

J = 1 

where 7 is the Euler constant and the notation is as bellow. 


1 . Introduction 


Let us consider the xi-function ^(s) = s{s — l)r(s/2)7r and the Li coefficients 

(A„)„>i defined by 

1 d” 


An — 




Then, the Li criterion says that the Riemann Hypothesis holds if and only if the coef¬ 
ficients (An) are positive numbers. Bombieri and Lagarias obtained an arithmetic 
expression for the Li coefficients An and gave an asymptotic formula as n —>■ oo. More 
recently, Maslanka m computed A„ for 1 < n < 3300 and empirically studied the 
growth behavior of the Li coefficients. Coffey mm studied the arithmetic formula and 
established a lower bounded for the Archimedean prime contribution by means of series 
rearrangements using the Euler-Maclaurin summation. In a generalization of the 
Li criterion for functions F in the Selberg class was given, and in m an explicit formula 
for the Li coefficients associated to F was established. 


The object of this paper is to derive a precise asymptotic formula for the generalized 
Li coefhcients using the saddle-point method. 
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The Selberg class S consists of Dirichlet series 


satisfying the following hypotheses. 


• Analytic continuation: there exists a non negative integer m such (s—l)™F(s) 
is an entire function of finite order. We denote by mp the smallest integer m 
which satisfies this condition. 

• Functional equation: for 1 < j < r, there are positive real numbers Qp, Xj 

and there are complex numbers fXj, to with > 0 and |w| = 1, such that 

(/)f(s) = - s) 


where 

r 

= F{s)Qp r(AjS + iij). 
i=i 

• Ramanujan hypothesis: a(n) = 0{rf). 

• Euler product: F{s) satisfies 


F(s) = ]^exp 

p 


jyks I 


with suitable coefficients b(p^) satisfying b(p^) = 0(p^^) for some 0 > g 
It is expected that for every function in the Selberg class the analogue of the Riemann 
hypothesis holds, i.e, that all non trivial (non-real) zeros lie on the critical line 5ft(s) = 
The degree of F G 5 is dehned by 


dp = 2^ Aj. 
j=i 


The degree is well defined (although the functional equation is not unique by the Le¬ 
gendre’s duplication formula). The logarithmic derivative of F(s) also has the Dirichlet 
series expression 

- —(s) = y] K(s) > 1, 

n—1 

where Ap{n) = b{n) logn is an analogue of the Von Mongoltd function A(n) defined by 


A(n) 


logp if n = p^ with fc > 1, 
0 otherwise. 


If Np{T) counts the number of zeros of F{s) G 5 in the rectangle 0 < 
|Q'(s)| < T (according to multiplicities), one can show by standard contour 
the formula 

NpiT) = ^TlogT + cpT + O(logr), 


3fi(s) < 1, 

integration 
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in analogy to the Riemann-Von Mangoldt formula for Riemann’s zeta-function C(s)! the 
prototype of an element in S. For more details concerning the Selberg class we refer to 
the survey of Kaczorowski and Perelli [6] . 


2. The Li criterion 


Let T be a function in the Selberg class non-vanishing at s = 1 and let us define 
the xi-function ^f(s) by ^f(s) = s™^(s — 1 )’”^(/)f(s). The function 5 f(s) satisfies the 
functional equation ^f(s) = w^f(1 ~ s). The function ^f is an entire function of order 
1. Therefore by the Hadamard product, it can be written as 

Cf(s) = Cf(0) n -) ’ 


where the product is over all zeros of Cf(s) in the order given by |5(/5)| < T for T —>■ oo. 
Let Af('r), n S Z, be a sequence of numbers defined by a sum over the non-trivial zeros 
of F{s) as 


AfH = 


1 -( 1 -- 
P 


where the sum over p is 


= hm 

^ T^oo ^ 
p |ap|<r 


These coefficients are expressible in terms of power-series coefficients of functions con¬ 
structed from the ^F-function. For n < — 1, the Li coefficients Af(r) correspond to the 
following Taylor expansion at the point s = 1 


dz 


+00 


log^F ) = Y 


n=0 


and for n > 1, they correspond to the Taylor expansion at s = 0 

|:logeF(^)=5:AH-»-iL”, 

^ n—0 


Let Z the multi-set of zeros of ^f{s) (counted with multiplicity). The multi-set Z is 
invariant under the map p i —> 1 ~ P, We have 



and this gives the symmetry Xpi—n) = Af(r). Using the corollary in O Theorem 1], 
we get the following generalization of the Li criterion for the Riemann hypothesis. 


Theorem 2.1. Let F(s) be a function in the Selberg class S non-vanishing at s = 1. 
All non-trivial zeros of F{s) lie in the line 5ie(s) = 1/2 if and only iflR. (Af(r)) > 0 for 

n = 1,2... 
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Next, we recall the following explicit formula for the coefficients XF{n). Let consider 
the following hypothesis: 

"H : there exist a constant c > 0 such that F(s) is non-vanishing in the region: 


•|s = (T-|-it; cr>l 


log(QF + 1 + \t\) 


Theorem 2.2. Let F(s) be a function in the Selberg class S satisfying LL. Then, we 
have 


Xpi-n) = mp+n{logQF - 


- E(n 


(- 1 ) 


1-1 




/7 

(/ - 1)! X—^+c 


E 

k<X 


^pik) TOF,, 

—^(logfc)' (logX)' 


r / +00 

+ ’^E^4"x7T77: + E 

4 = 1 

r n +00 

- 

j=l k=2 1=0 

where 7 is the Euler constant. 


+ Lj 


( 2 . 1 ) 


+ Lj + Lj) ^ 

1 ^ 

I F Xj F fJij 


Examples: 

- In the case of the Riemann zeta function, = 1, = 1, Ai = 1/2, and 

/ii = 0. With the equality 



(-i)‘ ( 1 -^) m. 


we hnd which was established by Bombieri and Lagarias Page 281]. 

- For the Hecke L-functions, Qp = nip = 0, Ai = 1, and fii = i, we find Af(?t-), 
which was established by X. -J. Li [9l p. 496]. 


3. Saddle-point method and the Norlund-Rice integrals 

Given a complex integral with a contour traversing simple saddle-point, the saddle- 
point corresponds locally to a maximum of the integrand along the path. It is then 
natural to expect that a small neighborhood of the saddle-point might provide the dom¬ 
inant contribution to the integral. The saddle-point method is applicable precisely when 
this is the case and when this dominant contribution can be estimated by means of lo¬ 
cal expansions. The method then constitutes the complex-analytic counterpart of the 
method of Laplace for evaluation of real integrals depending on a large parameter, and 
we can regard it as being 

Saddle-point method=Choice of contour -|- Laplace’s method. 

To estimate F{z)dz, it is convenient to set F{z) = where f{z) = fn{z), involves 
some large parameter n. We chose a contour C through a saddle-point 77 such that 
f'{r]) = 0. Next, we split the contour as C = and the following conditions are 
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to be verified. 

i) On the contour the tails integral/^(i) is negligible 

J F{z)dz = O (^J F{z)dz^ . 

ii) Along a quadratic expansion, 

fiz) = fiv) + ^f"{v){z - vf + 0{<pn) 

is valid, with (/)„—>■ 0 as n —>■ oo, uniformly with respect to z G 

ii) The incomplete Gaussien integral taken over the central range is asymptotically equiv¬ 
alent to a complete Gaussien integral with e = ±1 and /"(i?) = e*‘^|/"(? 7 )|: 


f + OO 


Jc(0) d-oo ' 

Assuming i), ii) and iii), one has, with e = ±1 and arg(/"(? 7 )) = (p 

o-iv/2 


2'K 


\nv)\' 


1 


27T 


V I/"WI 


=±, 


2'k 


!/"(??) I 


cfiv) 


This method is the main tool to prove our result. We finish this section by reviewing the 
definition of the Ndrlund-Rice integral. 

Lemma 3.1. Let f{s) be holomorphic in the half-plane 5ft(s) >170 — 5 . Then the finite 
differences of the sequence {f{k)) admit the integral representation 


E (fc)(-i)V(fc) = 

k—UQ 


i-iy 


fi 


2*71 Jc s{s — l)....{s — n) 


ds. 


where the contour of integration C encircles the integers {no, ..,n} in a positive direction 
and is contained in JR(s) > *70 — 5 - 

Proof. The integral on the right is the sum of its residues at s = no,.., n, which precisely 
equals the sum on the left. □ 

4. Asymptotic formula for the Li coefficients 

A natural question is to determine the asymptotic behavior of the numbers Xpin). 
Our main result in this paper is stated in the following theorem. 

Theorem 4.1. Let F{s) be a function in the Selberg class S. Then, under the Generalized 
Riemann Hypothesis, we have 


where 


Xpin) = -^n\ogn + cpn + O (Vnlogn) 


CF = ^(7 -1) + ^ log(Ag|), A = xf^ 

i=i 


and 7 is the Euler constant. 
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Remark 4.2. We conjecture that the asymptotic formula for the numbers Ap(n) in 
Theorem 14.11 holds for any function in the Selberg class without any assumption. 


(4.1) 


Four our purpose, it is sufficient to study sums of the form 

CG,f) 


H„(m,k) = J2(-iy(r)' 


1=2 




where m and k are tow complex numbers, (^(s, q) is the Hurwitz zeta function and defined 
for complex arguments s with 5ft(s) > 1 and q with ^{q) > 0 by 

+00 

Cis,q) = 


1 


n—0 


(n + g)®' 


Proposition 4.3. Hn{m,k), defined by (14.11) . satisfy the estimate 

Hnim,k) = ^ + log A: + 1 - hn-i^ + o„(m, fc), 

where the an{m,k) are exponentially small: 


a„(m, k) = 


_ 1 fY] 


1/4 


k \TrkJ 


Ann 


eM-\H-)cos W—-- - 


47 rn bn 2nm \ 




(r 


+ 0 


Here, /i„ = l + i + .. + i is a harmonic number, and ^/’(a;) is the logarithm derivative of 
the Gamma function. 

The previous results for the Riemann zeta function may be regained by setting m = 
k = 1 and for the Dirichlet L-functions m = 1 and k = 2. 

Proof. Convert the sum to the Ndrlund-Rice integral, and extending the contour to the 
half-circle at positive infinity. The half-circle does not contribute to the integral. One 
obtains 

(-1)” , C(s,f) 


Hn{m,k) = 


2in 


-ds. 


/ 3 / 2 -»oo A:'s(s - l)..(s - n) 

Moving the integral to the left, one encounters single pole at s = 0 and a pole at s = 1. 
The residue of the pole at s = 0 is 

k 


Res{s = Q) = C{{),-) = -—Y 


1 


kn 


sin 


1=1 


^nlm^ 


if 


= -Bi 


m 

T’ 


where ip is the digamma function, Bi is the Bernoulli polynomial of order 1, and 

71 / 777 - \ 

Res{s = 1) = - [ip{—) + \ogk + 1 - . 


Then we obtain 


Hn{m,k) = “ 7 -I- log/c -I-1 - hn-i^ + an{m,k), 

an{m,k) = O (e-"^) , 


where 
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for a constant K of order m/k. Indeed, we have 

^—1 /2 + ZCXD 


(-1)^ 

an{m,k) = n\ 


C( 5 ,f) 


'-1/2-ioo kh{s-l)..{s-n)' 


Recall that the Hurwitz zeta function satisfies the following functional equation 






k J [kirky 


cos 


1=1 


TTs 27rZm\ 

Y ~) 


C 


Therefore, 
an{m, k) 

(4.2) 


7Z\ V 
OhiiT 


r 3 / 2 -e*oo ^ r(s)r(s-l) /tts 2 T:lm\ f I 

■ cos I - — j C I s 


2 kiTr ^ ( 27 r)« r(s + n) 


ds 


Til K 
O h i'TT 


.3/2+,oo ^ r(s)r(s-i) / i 


2 kiiT 

n\ 

2 ki 7 T 


1^1 

k 

E 


/3/2-ioo ( 277 )'* r(s + n) 


C s, 


ds 


1 r(s)r(s-i) i 


/3/2-ioo ( 277 )* r(s + n) 


e* 2 C j s — ds. 


For large values of n, those integrals will be evaluated by means of the saddle-point 
method. Noting that the integrand in (14.21) have a minimum, on the real axis, near 

s = ao = \J and so the appropriate parameter is z = sjy/n. Changing s by z, and 
take z constant and n large. Then 

(4.3) ar!,(m,fc) = 

We have 


z=i 


^<70+200 


J{^)dz 


|•ao+ioo _ 
' aQ — ioo 


f{z) = logn! + -logn + (j){z^/n), 


with 

(t>{.s) = -slog 
using the approximation 

for larg s. Furthermore, 


2Trl 


TTS , 

- I— + log 


(^{s,iik) = {k/iy + 0 


r(g)r(s-i) 
F(s -I- n) 

I 


O 


k + l 


k + l 


+ 00 


logc(s) = E 


n —2 


Hn) 

n® log n ’ 


where A(n) is the Von-Mangoldt function. The asymptotic expansion for the Gamma 
function is given by the following Stirling expansion 

1 1 

logr(x) = {x- -)loga;-a:-k -log( 27 r) + ^ 


2 / 




2j(2j - l)x^i 


-1 > 
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where are the Bernoulli numbers. Expanding to 0(l/n) and collecting terms, we 
deduce 

' 2TrB 
~k 


1 


fiz) = Ti^ogn - zy/n [log 


t- + 2 - 21og2: ) + log(27r) 


1 


1 




(10 + 2") 

V 

^ 73 \ 


O 




2 6 ' 72zV 

The saddle-point is obtained by solving the equation f'(z) = 0, and we have 

zo = (1 + i)y 

Let a = a{n) = z^^ =(!-(- i)^^ which thus is also an approxi¬ 
mate saddle-point for Now, we fix the contour of integration and provide final 

approximations. The contour goes through the saddle point a = cr(n) and symmetrically 
through its complex conjugate a = a(n). In the upper half-plane, it traverses cr(n) along 
a line of steepest descent whose direction, as determined from the argument of 
The contour also includes parts of two vertical lines of respective abscissae Re{s) = a^/n 

and C 2 ^/n, where 0 < ci < < C 2 < 2^The choice of the abscissae, ci and 

^ ttZ > 


> 


C 2 , is not critical (it is even possible to adapt the analysis to ci = C 2 = y ^). 

• Contribution of the vertical parts of the contour are O , for some Lq 


2TTlr. 

k 


• The slanted part of the contour is such that all the estimates of /(z), /'(z) and 
/"(z) apply. The scale of the problem is dictated by the value of which 

of order We have 


/ =[ (l + 0 (n-^/^)). 

Jslanted ^central ^ ' 


/slanted ./central 
Using that zq = f{-/nzo) = —2^/nzo and 

f"{y/nzo) = 2/(Vnzo) -f 0{n~^). 
Now, we use the saddle-point formula to obtain 


./central 

Hence 

/ 

./slanted 


= g/(U" 2 o) 


/ 2. 

f 2T:^ln\ 

/ l/"(v^2o)| 

( k J 




/o f M I O , / 47r^Tr ^ 

T/Sgl, k ) 


O 
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Substituting, we obtain 
(4.4) 

r-o+^oo 


= f —-— j e*^ exp I —(1 -|- i) 


Anln 


+ O (n . 


The integral for / is the complex conjugate of (14.41) (having a saddle-point at the complex 
conjugate Zq). Finally, equations (14.3p and (14.41) together gives 


i{m,k) = 


1 (I 




1/4 


exp(- 


47rZn, 


cos 


47rZn Stt 27rZm \ 

k 8 k j 


+ O . 

For large n, only the I = 1 term contributes significantly, and so 


a„(m, k) = 


1 / 2 n\ 


1/4 


Ann. j Ann 5n 2nm\ / 

exp(-\/—)cos - - - —j + O ( 


k \nkJ 

which means that the terms a„ are exponentially small. 




□ 


Proof of Theorem 14.11 Since, we interest to the real part of Aip(n), then, without loss 
of generality, we assume that fij is a real number. First, write the arithmetic formula of 
Af(— n) (equation (12.If) ! as 

7 n 

Xri-n) = miP-Hn(logQiP-y7)-^(ni7E(^-l) 


Z =1 


(4.5) 

where 


- 

i=i V 


+ 00 

E' 


\j + fj,j 


k'j + A^j) y ^ 


(- 1 )' 

= “ir xEioo 


y/^(i„gt)<-iLL(i„g.Y)‘+' 


k<X 


l + l 


are the generalized Stieltjes constants and 


+ 00 




Note that 




+ 00 


^ = E(fc)(-^^)'E 


1 


(I Xj ,,v 

fc=2 1^0 ^ 3 ^ 3 

which, with the above notation of Hn{rn, k) (equation ()4.ip ). is equal to 

= il„(l + y,A-i). 


^ V ^ + Mi. 

Aj -l- fij) 
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Since Xj > 0 and fij = Re{fj,j) > 0, applying Proposition 14.31 with m = 1 + ^ and 


A: = we deduce 


(4.6) Ij — ^Xj + fij — —nXj {i>{Xj + ^ij) + log(Aj ^) + 1 — +a„(H- A^- ^), 


where 


f 2n _ / _ Stt a 

an(l + ^,A“^) = Aj AjJ exp(--\/47rnAj) cos |^y^47rnAj-^- 2Ti{Xj+^j)j 

+0 

The a„ are exponentially small, then 

(4.7) 

From (14.6|) and (14.7L we obtain 

(4.8) Ij = ^Aj + fij — —^ — nXj {'ipiXj + ^j) + log(A^ ^) + 1 — + 0(n). 

Summing (14.81) over j, we get 

(4.9) 

^ /j = ^ f Aj + /ij - - j - n'^Xj {ip{Xj + fij) + log{Xj ^) + 1 - hn-i} + 0{n). 
i=i i=i ^ ^ 1=1 

Using the expression 

+00 

1 Z 


an(l + ^>^j — C>(1)- 


ijj{z) = -1- - + J2 


z ^^Kl + z) 


where 7 is the Euler’s constant, and the estimate 

we deduce from (14.51) and (14.91) that 


Xpi-n) = ^Aj nlogn+< ^Aj (7-l) + log(5F + Xl^jlog^i 


K3 = i 


vi=i 


1=1 




Z =1 


Recalling that dp — Aj and noting that A = n X^^^ , we have 
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Now, we obtain a bound for SF{n) = — — 1) in terms of 

XF{-n,T):= ’ 


p; \Qp\<T ^ 

where T is a parameter. 

Lemma 4.4. If the Generalized Riemann Hypothesis holds for F G S, then 

Spin) = 0{\/n\ogn). 

Proof. The proof is an analogous to the argument used by Lagarias in We use a 
contour integral argument, and we introduce the kernel function 


fcn •— ( 1 + 


1 


-i = E(r) 


1=1 


If C is a contour enclosing the point s = 0 counterclockwise on a circle of small enough 
positive radius, the residue theorem gives 

I{n) = + = Spin). 

We deform the contour to the counterclockwise oriented rectangular contour C consisting 
of vertical lines with real part 5R(s) = uo and 5R(s) = ui, where we will choose —3 < uo < 
—2, (Ti = 2^JrL and the horizontal lines at S(s) = ±T, where we will choose T = y/n + en, 
for some 0 < £« < 1- The residue theorem gives 




— Sp{n) + E (1 + 


1 


P- 1 


-1 + 0 ( 1 ). 


p; |Sp|<T 

The term 0(1) evaluates the residues coming from the trivial zeros of F{s). Using the 
symmetry p \ — > 1 — p, we can write 


1 - 


- 1 = 


P- 1 


- 1 . 


Then 
We have 


I'{n) = Sp(n) - Xp(-n, T) + 0(1). 


\Xp{-n, \/n) - Xp(-n,T)l = O(logn). 

This follows from the observation that |T — ^/n\ < 1, that there are O(logn) zeros 
in an interval of length one at this height, and that for each zero p = j3 + i"f with 
\/n < l^(p)l < + 1 there holds 

V2 

< 2. 


fp-l\ 


1 

-— 

< 

1 + - 

\ p J 


n 
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We now choose the parameters (Tq and T appropriately to avoid poles of the integrand. 
We may choose tJo so that the contour avoids any trivial zero and T = ^/n + e„ with 
0 < Cra < 1 so that the horizontal lines do not approach closer than O(logn) to any zero 
of F{s). Recall from |16) that for —2 < 5R(s) < 2 there holds 

^(»)= E —+ o(iog(o,(i + i»i))). 


Then on the horizontal line in the interval — 2 < 5R(s) < 2, we have 

= 0 (log'T). 


F' 


The Euler product for F[s) converge absolutely for 5ft(s) > 1, hence the Dirichlet series 


for -^(s) converge absolutely for 5R(s) > 1. More precisely for a = 5R(s) > 1 


F' 

"f 

For a = 3?(s) > 2, we obtain the bound 


(cr) < OO. 


(a) < 2 -('"- 2 ). 


Consider the integral F {n) on the vertical segment (Li) having a\ = 2%Jn. We have 


F' 

— is) 

< 

F' 

F ^ ’ 


~F 


s 


< 


Then 


(l + ^) + 1 < (^1 + < exp(yH/2) < 2'/^. 

< C'o2"2(^+2). 


F' 


Furthermore, the length of the contour is we obtain | = 0(1). Let s = a + it 

be a point on one of the tow horizontal segment. We have T > so that 


1 + 


< 1 + 


cr + 1 


+ T 2 ’ 


By hypothesis > n, so for —2 < cr < 2, we have 


and 


|fcn(s)| < 1 


F' 


4 + n 


+ 1 = 0 ( 1 ) 


= 0(log2 T) = 0(log2n), 


since we have chosen the ordinate T to stay away from zeros of F{s). We step across the 
interval (+ 2 ) toward the right, in segments of length 1, starting from a = 2. Furthermore 


fcn(s + 1) + 1 


kn{s) + 1 


< 1 


1 


< e, 
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we obtain an upper bound for |^rt(s)-^(s)| that decreases geometrically at each step. 
After O(logn) steps it becomes 0(1), and the upper bound is 

=0(log^n + v^) =0(Vn). 

For the vertical segment (L 3 ) with 5R(s) = ctq, we have |fc„(s)| = 0(1) and |-|r(s)| = 
O [(5_F(log(|s| + 1))]. Since the segment (L 3 ) has length 0(i/n), we obtain 

\I'l^ = 0{y/n\ogn). 

Totalling the above bounds gives 

Spin) = A_f(— n, T) + 0{y/n\ogn), 


with T = y/n + e„. If the Generalized Riemann Hypothesis holds for F{s), then we 


have 


1 - 


p-i 


= 1. Since each zero contributes a term of absolute value at most 2 to 


Af(— n,T), we obtain using the zero density estimate {Np(T) ^ TlogT) 


Xp{-n,T) = 0{T\ogT + 1). 

Therefore Xp{—n, y/n) = 0{y/n\ogn), and Lemma lT4l follows. 


□ 


Using Lemma KM and the expression (I4.10p of Xp{—n) and Xp{—n) = Xp{n), we 
obtain 

^F(n) = ^nlogn+ 1 ^( 7 - 1) + ^ log(AQ|)| n + O(ynlogn), 

which concludes the proof of Theorem 14.11 
Exemples : 

- In the case of the Riemann zeta function, we have d(^ = 1, Qd; = and A = 1/2. 

This proves again under the Riemann Hypothesis the asymptotic formula established by 
A. Voros in [171 equation (17), p. 59]. 

- Also, in the case of the principal L-function L(s,7r) attached to an irreducible cuspidal 
unitary automorpohic representation of GL{N), as in Rudnick and Sarnak |141 §2], we 
have Dp = N, Qp = ( 3 ( 71 ) 71 “^/^, and A = 2“". We find under the Generalized Riemann 
Hypothesis the asymptotic formula for A„(7r) established by Lagarias in [T] equations 
(1.12) and (1.13), p. 4]. 
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